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ABSTRACT 

This paper aims to study the use of future SNAP data together with the 
result of searches for strong gravitational lenses in future large quasar surveys to 
constrain the Generalized Chaplygin Gas (GCG) model, with equation of state 
p = —A/p a , where A is a positive constant and < a < 1. The GCG is 
considered as a possible unification scheme for dark matter-dark energy. It is 
found that both experiments should be able to place important constraints on 
the model, especially when both tests are used together. 

Subject headings: cosmology: cosmological parameters, dark matter — equation 
of state — gravitational lensing 

1. Introduction 

In the last few years there has been mounting evidence that not only visible matter 
is a small component of the observable Universe, but that matter itself, dark or visible, 
makes up only a fraction of the cosmic energy budget. Indeed, SNe la experiments indicate 
that the Universe is expanding in an accelerated fashion (Riess et al. 1998; Garnavich et al. 
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1998; Perlmutter et al. 1999), and together with nucleosynthesis constraints (Buries, Nollet, 
& Turner 2001), and the CMBR power spectrum (Balbi et al. 2000; de Bernardis et al. 
2000; Jaffe et al. 2001) it is possible to arrive to a concordance where the Universe is made 
of about 5% of baryonic matter, 25% of dark matter and 70% of dark energy, a negative 
pressure component. These results are also in agreement with large scale structure (Peacock 
et al. 2001), and independent determinations of the matter density (Bachall & Fan 1998; 
Carlberg et al. 1998; Turner 2000). The nature of the dark matter and of the dark energy is 
not well established; dark matter candidates include axions, neutralinos and possibly a self- 
interacting scalar-field (see e.g. Bento, Bertolami, & Rosenfeld 2001 and references therein), 
while the most obvious candidate for the dark energy is the vacuum energy, or an uncanceled 
cosmological constant (Bento & Bertolami 1999; Bento, Bertolami, & Silva 2001). Another 
possibility is a dynamical vacuum (Bronstein 1933; Bertolami 1986a,b; Ozer & Taha 1987) or 
quintessence. Quintessence models most often involve a single scalar field (Ratra & Peebles 
1988a,b; Wetterich 1988; Caldwell, Dave, & Steinhardt 1998; Ferreira & Joyce 1998; Zlatev, 
Wang, & Steinhardt 1999; Binetruy 1999; Kim 1999; Uzan 1999; Amendola 1999; Albrecht 
& Skordis 2000; Bertolami & Martins 2000; Banerjee & Pavon 2001a,b; Sen & Sen 2001; 
Sen, Sen, & Sethi 2001) or two coupled fields (Fujii 2000; Masiero, Pietroni, & Rosati 2000; 
Bento, Bertolami, & Santos 2002). However, these models are afflicted with a fine tuning 
problem to explain the cosmic coincidence problem, that is, why did the dark energy start 
to dominate the cosmological evolution only fairly recently. 

Recently a new possibility has been suggested in Kamenshchik, Moschella, & Pasquier 
(2001), Bilic, Tupper, & Viollier (2002), and Bento, Bertolami, & Sen (2002a), which uses an 
exotic equation of state, and considers the evolution of the equation of state of the background 
fluid instead of a quintessence potential. This is achieved by using a background fluid, 
the generalized Chaplygin gas (GCG), described by the equation of state (Kamenshchik, 
Moschella, & Pasquier 2001; Bento et al. 2002a) 



The case a — lis called the Chaplygin equation of state, p c h = —A/p c h, while when a = 0, 
the model behaves as the flat ACDM, as discussed in Section II. The quantity A is related to 
the present day Chaplygin adiabatic sound speed through = aAj p c h,o, where p c h t o is the 
GCG density at present. This model allows for an unification of dark energy and dark matter 
and admits a brane interpretation (Bento et al. 2002a). Furthermore, analyses indicate that 
it may be accommodated within the standard structure formation scenario (Bilic et al. 2002; 
Bento et al. 2002a). These results have been challenged in Sandvik et al. (2002), however 
their analysis does not take into account the effect of baryons, which is expected and shown 
to be important (Bega et al. 2003), and makes the GCG model compatible with the 2DF 
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mass power spectrum. Also, Sandvik et al. (2002) analysis is based on the linear treatment 
of pertubations close to the present time, thus neglecting any non-linear effects. Thus we 
consider that a more careful and thorough analysis must be done before ruling out the GCG 
model on the grounds of its effects in the matter power spectrum or structure formation. 

Since the GCG exhibited great potentialities, it has been recently the subject of great 
interest, and various attempts have been made to constrain the model using the available 
observational data. In Avelino et al. (2002) the GCG was tested against the SNe la data and 
the matter power spectrum, and it is concluded that the GCG model agrees with the data, 
and interesting constraints on A s = Aj 'p^o are f° un d for fixed values of a. Also, the value 
A s — 1, corresponding to a pure cosmo logical constant, is strongly ruled out by data. The 
limits were somewhat degenerate mainly because three parameters, a, A s , and the density 
relative to critical of the Chaplygin gas, VL chl were considered. In Dev, Alcaniz, & Jain 
(2002) tests on the age of the Universe and strong lensing statistics were used to constrain 
the Chaplygin equation of state (a = 1), and it is found that the available data agrees with a 
Chaplygin gas model; from the age test it is concluded that A s > 0.96, and from the statistics 
of lensed quasars it is concluded that A s > 0.72. The age test together with SNe distance 
measurements were used to constrain the generalized Chaplygin equation of state in a flat 
Universe made up of the GCG and baryons, with fixed densities relative to critical, yielding 
confidence regions in the (A s ; a) plane (Makler, Oliveira & Waga 2002). The data indicated 
that 0.6 < A s < 0.85 at 95% confidence level (CL) and almost no constraints of interest 
to parameter a. The strongest constraints were obtained by Bento et al. (2002b) using the 
CMBR power spectrum measurements from BOOMERANG (de Bernardis et al. 2002) and 
Archeops (Benoit et al. 2002), together with the SNe la constraints from Makler et al. (2002). 
It is found that 0.74 < A s < 0.85, and a < 0.6, ruling out the pure Chaplygin gas model. 
Using the bound arising from the age of the APM 08279+5255 source, which is A s > 0.81 
(Alcaniz, Jain, & Dev 2002), fairly tight constraints are found, namely 0.81 < A s < 0.85, 
and 0.2 < a < 0.6, which also rules out the ACDM model. These results were found to 
agree with the WMAP data (Bento et al. 2003). 

The prospect of launching the SNAP satellite 1 will give rise to a new era of precision 
cosmology, and together with data from CMBR experiments and from large quasar and 
galaxy surveys, cosmologists will be able to test several models with unprecedented accuracy. 
Of course, this increase in accuracy may find its limits in new degeneracies, given the growth 
of free parameters, which in turn will require the use of more cosmological tests to distinguish 
competing models. In this paper we aim to study how well the GCG model can be constrained 
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with future data from the SNAP satellite, and from strong lensing statistics. We will be 
especially interested in the use of the Sloan Digital Sky Survey 2 (SDSS) data to constrain 
cosmological parameters, and also smaller surveys such as the 2dF 3 . Even though it has 
already been verified that the SNAP satellite will be able to constrain the model quite 
accurately (Avelino et al. 2002; Makler et al. 2002), the tighter constraints provided by 
the CMBR (Bento et al. 2002b, 2003) invite us to go a step forward as SNe confidence 
regions depend strongly on the fiducial model considered. Moreover, we add that the use 
of gravitational lensing statistics from future large surveys has not been explored in the 
context of the Chaplygin gas. Thus, the main aim of the our study is to use all three tests 
to constrain as much as possible the parameters of the GCG. 

The outline of this paper is as follows. In Section 2 we explain how the Chaplygin model 
closely mimics the ACDM model, and show that the case a = corresponds to the ACDM 
model. In Section 3 we consider the magnitude versus redshift cosmological test. In 3.1 we 
briefly describe the nature of this test, and in 3.2 we estimate the errors considered. In 3.3 
we describe the method used to obtain the confidence regions shown in 3.4. 

Section 4 is dedicated to strong lensing statistics. In 4.1 we briefly describe the strong 
lensing test, in 4.3 we explain how we use the Fisher matrix formalism to build the confidence 
regions that are exhibited in 4.4. 

2. Generalized Chaplygin gas compared to flat ACDM. 

The cosmological tests considered here depend on the cosmological model through the 
comoving distance r(z, H ,&), where is a vector containing the cosmological parameters 
we are interested in determining, and H Q being the Hubble constant. For a Friedman- 
Robert son- Walker (FRW) Universe the comoving distance as function of redshift is given 



by 





2 www. sdss.org 
3 www. 2dfquasar.org 
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where Q k = l — Q being the total density of the Universe relative to the critical one. Only 
flat cosmological models, f4 = 0, are considered in this paper. 

The expansion rate may be written in terms of the total mass-energy density as a 
function of redshift. Recalling the Friedman equation, 

H\z) = ^p{z) , (3) 
where p(z) is the total mass-energy density at redshift z, the expansion rate is given by 

where po = p(0) is the energy density at present. For a flat Universe, po = p c , where p c is 
the critical density. 

For the flat AC DM, ignoring the small radiation energy density contribution, one has 
p(z) = po [Q A + (n CD M + to b ) (1 + zf] , (5) 

where Qqdm, are the baryonic, cold dark matter and vacuum energy densities relative 
to critical, respectively. Thus, the Hubble parameter becomes 

H(z) = H [n A + (Q CDM + Q b )(l + zf] 1/2 . (6) 

The GCG unifies dark matter and energy, but does not account for the baryons. Therefore in 
a consistent model of the Universe baryons and radiation must be included. The contribution 
of radiation shall be ignored, since its effect is not important for the redshifts considered 
here. Hence, for a Universe made up of the GCG plus baryons one has (Bento et al. 2002a) 



Pch{Z) 



p 6 ,o(l + zf +(A + B(l + z f^)) 1/{1+a) . ( 7 ) 



where p^o is the baryon density at present, and B is an integration constant. Using the 
freedom to define the integration constant B, this expression may be rewritten as 

p ch (z) = p 6 , (l + zf + p cM [A. + (1 - A s )(l + zf 1+ ^] 1/(1+a) , (8) 

where p c hfl = (A + is the present day density of the Chaplygin gas, and A s = 

A/pIho- Thus, for a flat Universe, the Friedman equation becomes 

H ch (z) = H {Q b (l + zf + Q ch [A s + (1 - A s )(l + ^)3(i^)] 1 /d+-)| 1 /2 ^ ^ 
where Q c h = p c h,o/Pc is the GCG density relative to the critical one. 
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Even though the generalized Chaplygin equation of state looks rather exotic, the result- 
ing cosmological model behaves in a manner that mimics the ACDM model (Bento et al. 
2002a). First let us consider the very early Universe, when its size was very small compared 
to the present one, a(z) <C a = 1, or z ^> 0. The ACDM behaves as a flat CDM dominated 
Universe: 

p( z > 0) = po (ficoM + n b ) (1 + zf . (10) 

The GCG also behaves as a flat CDM dominated Universe with a dark matter density 
&cdm = ^ch(l ~~ A s y/( 1+a \ as may be seen in 

p ch (z » 0) = po (Q b + Q ch (l - A s ) 1/(1+q) ) (1 + zf . (11) 

This property ensures that the GCG model is consistent, at large scales, with the CDM 
structure formation scenario. 

Consider now the late Universe, when a ^> a . The ACDM behaves as a vacuum 
dominated Universe, 

p{a > a ) = Po^a , (12) 

while the GCG model also behaves as a vacuum dominated Universe, that at present has a 
vacuum density = Q c hAl^ 1+a \ 

p ch (a » ao) = po^cU s V(1+Q) • (13) 



Between these two limits, the GCG density may be expanded in subleading terms, 
yielding 



p ch = po < Q b (l + zf + tt ch 



A]l^ + - 1 1 (1 + z f 1 ^ 
1 + a A a/{1+a) 



(14) 



Thus, between the dust and the De Sitter phases, the GCG behaves like a two component 
fluid made up of vacuum energy with a present day density given by Q\ = Q c hAl^ 1+a \ and 
soft matter with the equation of state p = ap, and at present 

1 \ 1-As 



a 



ch 



1 + aJ A 



<a/(l+a) 



(15) 



Also, making the associations 



^ a A s 



nira-As) 



1+Q 

CDM > 



(16) 
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and using them in Eqs. (8-9), one obtains 

Pch(z) = p fe , (l + ^) 3 + Po^r + ^^ 

H ch (z) = H {sk{l + + K +a + &(1 + ^) 3(1+a) ] V(1+Q) } 

[n\ + " + n}% M ] 1/{1+a) + n b = i, 

which allows one to see that a = is identical to the flat ACDM model. 



1/2 



(17) 

(18) 
(19) 



3. Magnitude versus redshift cosmological test. 

3.1. Description of the test. 

The magnitude versus redshift test explores the dependence of the distance modulus of 
a source on its redshift. For a source with redshift z, its apparent magnitude is related to 
its absolute magnitude through 



m(z) = M + 5\ogd L (z,H ,@) + 25 



(20) 



where represents the several cosmological parameters we are interested in, di is the lumi- 
nosity distance, d^ — (1 + z)r, measured in units of 10 pc, and r is the comoving distance, 
Eq. (2). The luminosity distance is then given by 



d L (z,H ,Q) 



(l + z)c 
Ho^—Qk 



:i+z)c 



sin 



H^—Vt k dz' 
o H(z',H ,Q) 

dz' 



o H(z',H ,Q) 
sinh 



(l + *)c _,_. n . ( r H Q ^/rr k dz' 



o H(z',H ,Q) 



tt k < , 

fi fc = o, 
n k > o . 



(21) 



Defining a new quantity D L = H d L one may group the terms that depend on the 
absolute magnitude and H . The apparent magnitude then reads 



m(z,H ,e) = M-51og(# ) + 25 - 5\ogD L (z,Q) 
= M-5\ogD L (z,Q) . 



(22) 



The quantity M. = M — 5 log H + 25 is usually called zero point value or intercept since it 
it the value of the apparent magnitude at the point where \ogD L = 0. This quantity can be 



- 8- 



measured using low redshift SNe, or be treated as a statistical nuisance that is marginalized 
(Permutter et al. 1997). We shall follow Goliath et al. (2001) and assume either an exact 
knowledge of M. or no prior knowledge of it. This allows us to concentrate on the study of 
the parameters that are relevant to us. 



3.2. Error estimates for one year of SNAP data. 

In the next subsections we aim to study the ability of the SNAP mission to test the GCG 
equation of state. We shall use the expected error estimates for the SNAP satellite from 
Weller & Albrecht (2002), and consider that the systematic errors for the aparent magnitude 
m are given by 

/ x 0.02 

<J S ys{z) = ' ( 23 ) 

which are measured in magnitudes such that at z = 1.5 the systematic error is 0.02 mag- 
nitudes, while the statistical errors for m are estimated to be a s t a = 0.15 magnitudes. We 
place the SNe in bins of width Az pa 0.05. We then sum both kind of errors quadratically 



<Tma B (z i ) = J< S (^) + ^f > ( 24 ) 
y <n 

where rii is the number of supernovae in the redshift bin. The distribution of supernovae in 
each redshift bin is the same considered in Weller & Albrecht (2002), and is shown in Table 
1. 



3.3. Confidence regions for SNe la tests. 

To build the confidence regions a fiducial model corresponding to a vector S trU e is cho- 
sen, and log-likelihood functions x 2 are calculated based on hypothetical magnitude mea- 
surements at the various redshifts. This function is given by 

X 2 (M 0) = ]T H^'M ) -m(z,M,® theory)f ^ 
zi=0 ° ^ 

where the sum is made over all redshift bins and m(z,A4, 0) is defined as in Eq. (22). 

The x 2 function is directly related to the maximum likelihood estimator of a normal 
random variable, L = exp(— x 2 /2). Suppose one wants to impose a Gaussian prior to one of 
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the parameters being measured, 9, centered around 9o, with variance aj. Using the Bayes 
theorem one finds that 

1 2 , f(0-0o) 2 



^ = exp(-- X 2 )exp(^-^) . (26) 



Thus, the y 2 function is changed to 

2 2 , (0 ~~ ^o) 2 / -n 

X -»• X + 2 • ( 2T ) 

If one aims to impose a prior on a given parameter one is not measuring, one has to 
integrate the parameter out of the likelihood function using its probability density function, 
p(9), as prior, that is 

4rior on e = ~ 2 ln M ex P ("^) p(*) • (28) 

A special case is the marginalization over the zero point M.. When this zero point is 
exactly known it cancels out from the \ 2 expression, and the modified log-likelihood function 
X 2 (Goliath et al. 2001) is obtained, 

n A 2 
i=i 1 

where 

A = 5 \og w [D L (z, &)} - 5 \og w [D L (z, & true )\ . (30) 

Assuming no prior at all about Ai, one has to integrate it out of the \ 2 function, 
obtaining a modified log-likelihood function x 2 '- 

~ f +co ( 1 

X 2 = -21n / dMexp[--x 2 



— oo 
2 



r 



C \2tc 

where 



= r-^r • In(^) , (31) 



B = J2^, (32) 



C = Ei- (33) 
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3.4. Expected Confidence regions for 1 year of SNAP data. 

The values used for each fiducial model are shown in Table 2. As mentioned above, 
the use of SNAP data to constrain the GCG model was already considered in Makler et al. 
(2002) and Avelino et al. (2002), but our analysis differs these ones in three points. Firstly, 
the Avelino et al. (2002) analysis is done in the (0, ch: A s ) plane for fixed values of o, while 
we consider fl c h fixed, and find constraints in the (A s ,a) plane. Secondly, in Makler et al. 
(2002) only a ACDM fiducial model is considered, essentially with the purpose of testing 
whether SNAP will be able to rule out the GCG model, noting that it will be possible to 
distinguish between a ACDM and the Chaplygin gas (a = 1) model. Our analysis not only 
considers the capability of ruling out the GCG model, via our study of Model I, but also 
considers the precision of the constrains one may find for both our fiducial models. The 
inclusion of this second fiducial model is important, as it corresponds to the center point 
of the range allowed by CMBR (Bento et al. 2002b) and the age test (Alcaniz et al. 2002). 
Thirdly, we are not only interested in the use of SNAP to constrain the model, but also the 
use of strong gravitational lensing statistics, and especially to study what constraints might 
arise when both probes are used together. Throughout the paper we will also use the CMBR 
constraints from Bento et al. (2002b), essentially with two purposes. Firstly to compare 
these available constraints with the future ones, to see whether future tests will be able to 
improve the available constrains or merely confirm them, and secondly as a third test that 
might be used to improve the confidence regions. 

The expected confidence regions that are obtained are shown in Figures 1-3. In all cases 
a fixed baryon density with value = 0.05 is used. Two fiducial models were considered. 
Model I corresponds to the center value of the parameter range found by Bento et al. (2002b), 
(A 3 , a) = (0.83, 0.4). Model II corresponds to a flat ACDM model with Q A = 0.72, Q C dm = 
0.23, and = 0.05. Using the associations given in section 3.1, this corresponds to (A s , a) = 
(0.758, 0). Both sets of values are compatible with the CMBR data (Bento et al. 2002b) and 
available SNIa data (Makler et al. 2002). The set corresponding to the ACDM fiducial 
model is ruled out by the age estimate of the APM 08279-5255 source (Alcaniz et al. 2002), 
but it was included to illustrate the effect of choosing a fiducial model. 

The first conclusion one draws is that in order to achieve some precision, a good estimate 
of the intercept M is required. Without imposing prior knowledge of M. it is possible to 
distinguish only between a ACDM model (a = 0) and a generalized Chaplygin model with 
a > 0.75 at 95% confidence level (CL), or with a > 0.4 at 68% CL, as shown in Figures 2 
and 3. For Model I (Figure 2) the situation is even worse, and without a good knowledge of 
the intercept it is only possible to distinguish between models with a > 0.9 at 68% CL. In 
either case the pure Chaplygin gas (a = 1) is ruled out with 68% CL. To distinguish between 
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both fiducial models a good estimate of M is necessary. 

In the previous paragraph the ability of the test to distinguish between models was 
considered. Next we consider whether it is possible to improve on the limits imposed by the 
CMBR (Bento et al. 2002b), available SNe (Makler et al. 2002) and the age test (Alcaniz et 
al. 2002). We find that the answer depends on the fiducial model considered. For Model II it 
is possible to confirm or to reject the present limits. Indeed, using CMBR constraints from 
Bento et al. (2002b), Figure 1, leads to a smaller parameter region, but it allows no more 
than to assert that 0.2 < a < 0.65 and 0.81 < A s < 0.85 at 68% CL, or 0.15 < a < 0.75 
and 0.79 < A s < 0.87 at 95% CL. 

For Model I, the situation is different. This model is not compatible with the age 
estimate of the APM 08279-5255 source (Alcaniz et al. 2002), therefore the comparison 
must be made with the allowed parameter range from available SNe data and the CMBR 
constraints from Bento et al. (2002b). From only these two tests one can establish that 
0-74 < A s < 0.85, and < a < 0.6 (Bento et al. 2002b; Makler et al. 2002). For this fiducial 
model, the SNAP data will allow for tighter constraints, yielding 0.75 < A s < 0.80 and 
a < 0.35 at 95% CL. The use of additional CMBR constraints from Bento et al. (2002b) will 
allow further constraining a to a < 0.3 at 95% CL. 



4. Strong lensing probability. 

4.1. Description of the test. 

We start the description of the strong lensing probability test by setting the assumptions 
which are most commonly used in gravitational lensing statistics, namely: 

• Lensing galaxies are described by singular isothermal sphere (SIS) profiles. 

• There is no evolution in the population of lensing galaxies. 

• The Tully-Fisher and Faber- Jackson relations are independent of redshift. 

Under these assumptions, in a flat FRW Universe the probability that a source at redshift 
z s will be lensed by a galaxy is given by (Turner, Ostriker, & Gott 1984; Gott, Park, & Lee 
1989) 

r = F r -^- , (34) 
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where R = c/H , r(z) is the comoving distance, Eq. (2), and F is the dimensionless 
parameter 

f = 7h! < n ° a " > ' (35) 

no being the number density of galaxies at present and cry is their velocity dispersion. The 
average is computed over all lensing galaxies using the Tully-Fisher and Faber- Jackson re- 
lations together with the galaxy luminosity function (Fukugita & Turner 1991). We shall 
use the value estimated in Cooray (1999), F = 0.026. The information concerning the 
cosmological model enters through the comoving distance. 



4.2. Magnification bias and selection effects. 

To sucessfully use strong lensing statistics as a tool to test cosmological models, a 
detailed account of possible systematic effects is required. Here we provide a brief discussion 
of some effects and how to quantify them. 

One of the most important systematic effect is the magnification bias (Turner et al. 
1984). It affects magnitude limited surveys, since gravitationally lensed sources are pref- 
erentially included in the survey. This bias may be quite large, depending on the quasar 
luminosity function. Other important biases are the selection effects due to the limitations 
on the dynamical range (magnitude differences between the lensed images), resolution limi- 
tations and the presence of confusing sources such as stars (Kochanek 1991, 1993a). Besides 
these there is the issue of selection effects associated with quasar surveys (which tends to 
misrepresent the quasar population) such as the luminosity of the lensing galaxy (which 
tends to overwhelm the quasar luminosity and exclude it from the catalogue) and the red- 
dening due to the lensing galaxy (which may reduce the average magnification bias, increase 
flux ratios and absorb single images) (Kochanek 1991). 

Effects due to the lensing galaxy luminosity are expected to be marginal for bright 
quasar surveys, but for deep (low luminosity cut-off) surveys such as the SDSS the effect 
might be statistically dominant. Following Cooray & Hutterer (1999), we shall not include 
this correction. A correct estimate of these effects requires an accurate characterization of the 
apparent magnitudes of lensing galaxies as function of redshift (Kochanek 1991). Since there 
still is no consensus on how to quantify the reddening effect (Malhotra, Rhoads, & Turner 
1997; Falco, Kochanek, & Munoz 1998), we shall follow, as before, Cooray & Hutterer (1999) 
and disregard this effect too. 

The effect of confusing sources should be negligible if the lensing survey magnitude 
range is not excessively large (Am > 2.5 for a survey with M Um = 20 mag) (Kochanek 
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1993a). Thus, we shall disregard corrections of this nature, but remark that the necessity of 
including this correction depends on the survey characteristics. 

All these effects are somewhat connected, especially the magnitude bias, angular sep- 
aration and dynamic range effect (Kochanek 1993a,b), but since we are considering a hy- 
pothetical survey, we do not need to consider a detailed selection function. Hence we will 
follow the simplified approach of Cooray (1999), and write the probability that a source at 
redshift z will be lensed and detected as 

p(z) = tB(< m)f(A0) , (36) 

where B(< m) is the average magnification bias for the survey, which will be defined, and 
f(A9) is the angular selection function. Since we found no information on the expected 
magnitude difference limit for the SDSS we choose just to consider the correction due to 
the angular resolution. The ignored correction should be small (~ 0.9997 as estimated in 
Cooray (1999) for Am < 2.5, see also Kochanek (1993a)), hence ignoring this effect will not 
yield a large error. 

Considering the assumptions commonly used in the literature (see, for instance, Fukugita 
& Turner (1991)), the magnitude bias at a given magnitude level, B(m), is given by 

B{m) = N^mj r Nq{m + A ) P ( A ) rfA > ( 37 ) 

where Nq{m) is the intrinsic quasar number counts, A = 2.5 log A, A being the total mag- 
nification produced by the lens, and P(A)dA is the probability distribution of A. The 
probability distribution of the magnification A is given by (Turner et al. 1984) 

P(A)dA = -^dA A>2, (38) 

or alternatively 

P(A)dA = 7.37 x l(T a8A dA A > 0.75 . (39) 

We shall use the quasar number counts from Hartwick & Schade (1990), 
N (m ) - J (10/D°)10°- 86 (™- 19 - 15 ) m B < 19.15 , 

Q[ b) ~ \ (l /nO)l 0-28(m B -19.15) mR > 10 15 ^> 

The D represents the number of sources per square degree of sky. Note that the normal- 
ization of the number counts does not influence the magnification bias. The magnification 
bias is given by (Fukugita & Turner 1991) 

59.5 x io°- 06 ( 19 - 15 - m s) _ 59.2 mB < 18.40 , 

B(m B ) = { 2.5 x io°- 58 ( 19 - 15 - m s) 18 40 < m B < 19.15 , (41) 

2.50 m B > 19.15 . 



Since we are working with a hypothetical data set, we must average B(m) over the 
observed magnitude distribution. The average bias for a survey whose magnitude limit is 
ttil, B(< ttil), is then 



B(< m B ) 



f^N Q (m B )B(m B )dm B 
jT^ N Q (m B )dm B 



(42) 



For the SDSS, Cooray & Hutterer (1999) used a limiting magnitude of m L = 21 mag, 
although the survey has a limiting magnitude of m L = 23 mag. We shall follow this sugges- 
tion and also use = 21 mag, which corresponds to an average magnification bias factor 
of B(< m L ) = 2.89. 

Moreover we will use a step angular selection function, such that the lens is detected 
if the angular separation is between 0.1" and 6" , being undetected otherwise. For a SIS 
model, the fraction of lenses in this range is 0.901 (using the results from Kochanek (1993b); 
Fukugita & Turner (1991)). We shall consider a simplified selection function such that all 
lenses within this angular separation range will be recovered and write /(A#) = 0.901. The 
angular limits used are very optimistic, but this provides an estimate of how precise the 
survey must be to obtain sensible results. 

Clearly, what is observed in strong gravitational lensing statistics is the fraction of 
lensed sources relative to the overall number of sources. Consider that Nfotai sources (such 
as quasars or galaxies) are found. To study the effect that the size of the survey will have on 
the precision of the test, we consider three test models, with N Tota i = 25000, 10 5 , 10 6 sources 
(the SDSS survey is expected to yield the spectra of 10 5 quasars, and about 10 6 optical 
quasars. The 2dF survey expects to find 25000 QSO), spread uniformly in 40 redshift bins 
of width Az = 0.1, starting from z min = until z max = 4. In this model there will be 
n = NTotai/^O sources in each redshift bin. The effect of considering a coarser angular 
resolution may be seen as having less sources at each redshift. Thus, the 10 6 case may be 
seen as the best scenario. A more realistic angular selection function would give results that 
are between the 10 5 and the 10 6 sources cases. 

Consider now a fixed redshift bin Z{. For this redshift bin n = A^ oto ;/40 sources are 
picked, each one with an independent probability of being lensed given by p(zi). The number 
of lensed sources with redshift Zi, Ni(zi), is then a binomial random variable with parameters 
n and p(zi). Since n is large, and p(zi) is very small, one may use the Poisson PDF with 
average 



N exp (zi) = np(zi) . 



(43) 
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4.3. Building Confidence Regions. 

In order to build up confidence regions we employ the Fisher information matrix method 
(see e.g. Tegmark, Taylor, & Heavens (1997)), which provides an estimate of the inverse 
covariance matrix when the likelihood function of the random variable is known. 



The Fisher matrix F is given by: 



/d 2 \nL\ 

F * = - <75f7*r> ( 44 ) 



dOidOj 

where L is the likelihood of observing the data set x given the parameters 9±, . . . , Q v , and the 
average is being computed at the point corresponding to the true set of parameters @tme- 
The Fisher information matrix is the expectation value of the inverse of the covariance matrix 
at the maximum likelihood point, and it gives us a measure of how fast (on average) the 
likelihood function falls off around the maximum likelihood point. Assuming the Poisson 
statistics as a good description of the expected number of lenses as function of redshift, the 
likelihood function is written as 

L[n L (Az 1 ),n L (Az 2 ), n L (Az max )] = 

Zrnax -N exp (Azi) lU«l(A2i) 

= S MMV ■ (45) 

Taking the logarithm and calculating the derivatives one obtains (Cooray & Hutterer 1999) 

F \- 1 dN exp (z, Az) dN exp (z, Az) 
13 ^N exp (z,Az) 09 t 09, ■ 1 } 



Confidence regions are built by imposing a confidence limit on the likelihood function. 
To build these regions using the Fisher matrix one starts by expanding the logarithm of the 
likelihood in a Taylor series around the maximum likelihood point. Since this corresponds 
to a maximum of L, the first derivative vanishes and one finds 

2\n(L max /L)=6OF50, (47) 

where 56 = (6>i - 9 hmax ; ...;9 U - 0„, mox ). By choosing a range of values of 2ln(L max /L) 
corresponding to the desired confidence level, a second order equation in the 86i is obtained 
which can be solved to obtain the desired confidence regions. It can be shown that in the 
large n limit, 2 \n(L max /L) has a x 2 distribution with v degrees of freedom (Kochanek 1993b). 
For our study v — 2, and the 68.3%CL are found by imposing 2\n(L rnax / ' L) = 2.3, while the 
95%CL regions are found by imposing 2\n(L max /L) = 5.99. 
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It should be pointed out that this population model is not very realistic since fewer 
sources are found at low redshift due to the smaller volume involved and, on the other hand, 
fewer sources at high redshift as only very bright sources are visible. Considering a more 
realistic distribution should not change the results appreciably though, as one finds more 
sources at the medium redshift range, thus with a higher lensing rate, but on the other, fewer 
sources are found at high redshift, which results in an overall lower lensing rate. These two 
competing effects should thus cancel each other, and produce a small net effect. 

4.4. Results. 

The confidence regions that were built are shown in Figures 4-8. Since their implications 
are quite different, the results for each fiducial model are analysed separately. We consider 
strong lensing statistics as probe by themselves, and together with CMBR peak location and 
SNAP constraints. This goes beyond the analysis of Makler et al. (2002) where only the 
SNAP data was considered, and the analysis of Bento et al. (2002b, 2003), where only the 
CMBR peak location was analised. 

44.1. Model I. 

The gravitational lensing statistics test by itself is unable to yield important constraints 
on the parameter space of the GCG model, and, at most allows asserting that 0.76 < A s < 
0.94, with no limits on a. Even considering the CMBR peak position data from Archeops 
and BOOMERANG only an upper value a < 0.7, and 0.76 < A s < 0.87 is found, as can be 
seen in Figure 4. 

Thus, combining CMBR constraints from Bento et al. (2002b), lensing statistics and 
SNAP data (Figure 5) does not improve the limits found in Section 3.4. There is only a 
marginal improvement in the upper range of the allowed values of a. Therefore it is expected 
that, if the limits imposed in Bento et al. (2002b) are accurate, SNe distance measurements 
and strong gravitational lensing statistics we will not be able to substantially improve the 
present constraints, at least in the near future. 

44.2. Model II. 

For the second fiducial model, the results are more promising. As may be seen in Figure 
6, the search for lenses in the SDSS will allow to impose a fairly tight parameter range 
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a < 0.3, and 0.75 < A s < 0.85. For smaller lens surveys, these results together with CMBR 
peak location data may yield some constraints. It is particularly interesting to note that a 
survey with the size of the 2dF survey might be able to place some interesting upper limits 
on A s and a in the very near future, although it should be realized that using a more careful 
treatment of the selection effects might broaden the confidence regions. 

The most interesting results appear when we use both SNAP data and lensing statistics 
for a large survey such as the SDSS (Figures 7 and 8). We are close to determining the 
model, as they will allow to constrain the parameters to be in the range a < 0.1 and 
0-75 < A s < 0.77 at 68% CL, or a < 0.2 and 0.75 < A s < 0.79 at 95% CL. This means that 
SNAP data together with strong lensing statistics from a large precise survey should be able 
to rule out the CGC model as an unification scheme between dark matter and dark energy. 

4.5. Conclusions and Outlook. 

In this paper we have studied how the SNAP data and the data from future quasar sur- 
veys might be used to constrain a generalized Chaplygin equation of state. Our conclusions 
depend on the models under consideration. 

If the fiducial model is the model favored by available SNe (Makler et al. 2002), CMBR 
peak location (Bento et al. 2002b) and age estimates of distance sources (Alcaniz et al. 
2002) (Model I in Table 2), then the SNAP data will be able to definitely rule out both 
the pure Chaplygin well as the ACDM model, but it will not be able to provide 

tighter constraints than those already available. This is particularly so without a precise 
determination of the Hubble constant, H , and of the SNe absolute magnitudes, M in Eq. 
(22). In fact, to completely rule out the ACDM model, a good determination of these 
quantities is absolutely required. 

Using the available data on the CMBR peaks positions from BOOMERANG (de Bernardis 
et al. 2002) and from Archeops (Benoit et al. 2002), a slight improvement in the allowed 
parameter range is expected, but most likely to verify the range of parameters determined 
in Bento et al. (2002b), namely 0.81 < A s < 0.85 and 0.2 < a < 0.65 at 68% CL, or 
0.15 < a < 0.75 and 0.79 < A s < 0.87 at 95% CL. Of course that the situation will improve 
greatly with data from WMAP, especially with a more precise determination of the location 
of the acoustic peaks (Bento et al. 2002b). 

Lensing statistics by themselves will not be very helpful in constraining the allowed 
range of parameters, but together with SNAP data and CMBR constraints, the allowed 
parameter range is slightly tighter, yielding 0.81 < A s < 0.85, and 0.2 < a < 0.6 at 68% CL. 
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We therefore conclude that in order to achieve more precise limits on our fiducial Model I, 
the best strategy is to aim for better determinations on the location of the CMBR acoustic 
peaks. If the Universe is best described by Model I, the data from SNAP or from the search 
for lenses in large quasar surveys will only be able to verify the present day limits, and only 
marginally to improve them. 

Turning now to the AC DM fiducial model (Model II in Table 2), the results are much 
better. Again, an accurate determination of the zero point M. is required in order to achieve 
some precision. With such a determination one might rule out models with a > 0.2 at 68% 
CL, or with a > 0.35 at 95% CL. Using the CMBR peaks locations, marginal improvements 
may be obtained. The best results arise from the use of SNe and lensing statistics, where 
one comes close to determining the model, obtaining the allowed range, a < 0.1 and 0.75 < 
A s < 0.77 at 68% CL, and a < 0.2 and 0.75 < A s < 0.79 at 95% CL. The lensing results 
should be regarded with care, since our best case scenario is overtly optimistic. 

Therefore, if one considers a two component Universe made up of baryons and the GCG, 
future data from SNAP and from large QSO surveys, such as the SDSS, will be able to test 
the GCG model, and rule out a possible unification scheme between dark energy and dark 
matter it provides. The question of how well we will be able to constrain the parameter 
space depends on the model that is ruled out. 

4-5.1. Assumptions made for strong lensing statistics. 

As a final note, we mention that caution should be exerted when considering results 
from gravitational lensing statistics, as there are several systematic effects that we have not 
explored. The inclusion of a core radius on the galaxy model, for instance, will greatly 
suppress the lensing optical depth (Hinshaw & Kraus 1987), but on the other hand it will 
enhance the magnification bias (Kochanek 1996). This partly compensates the lower lensing 
probability and produces a small net effect. Also, the lack of central images in observed 
gravitational lenses indicates that the core radii will be smaller than 200 pc (Wallington 
& Narayan 1993). These results also agree with Kochanek (1996), and Fukugita & Turner 
(1991) where strong lensing statistics were used to constrain several systematic uncertainties. 
Therefore, the error from ignoring the core radius appears to be negligible. However, since 
one is dealing with an unprecedented degree of precision, uncertainties will inevitably affect 
the results. 

Effects due to galactic evolution have been studied in Mao (1991), Mao & Kochanek 
(1994), and Rix et al. (1994), with the conclusion that evolution will only affect the statistics 
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substantially if it is important at z < 1. For a vacuum dominated Universe this should 
not be a problem since galaxies form earlier that for a matter dominated Universe (Carrol, 
Press, & Turner 1992). Still, evolutionary studies of galaxies have not reached conclusive 
results. For instance, results from the Canada-France Survey (Lilly et al. 1995) indicate 
that the comoving number density of sources is independent of redshift. On the other hand, 
results from the Hubble Deep Field depend on the methods used to estimate the redshift 
distribution of sources. In Mobasher et al. (1996) it is argued that there are no signs of 
evolution, while Sawicki, Lin, & Yee (1997) indicate that there is some evolution, although 
not a dramatic one, from 0.2 < z < 0.5 until 1 < z < 2. On the other hand, a more dramatic 
evolution from z = 2.2 until z = 0.6 is argued in Gwyn & Hartwick (1996). 

In calculating the parameter F in Eq. (34), we must relate the velocity dispersion of 
matter in the galaxy to its luminosity. The problem of estimating the dark matter velocity 
dispersion ctdm from the luminous matter in elliptical galaxies arises, and is of great im- 
portance in gravitational lensing statistics, since a 22% change of cry in Eq. (34) leads to 
a 125% change in the expected number of lenses and a 50% change in image separations. 
Here we assumed that the velocity dispersion of dark matter is the same as the luminous 
matter, as considered in Kochanek (1991, 1993a,b, 1996), Chiba & Yoshii (1999), Cooray 
(1999), Cooray & Hutterer (1999), and Cooray, Quashnock, & Miller (1999). This equality 
was advocated in Kochanek (1993b) and Kochanek (1996). Other studies have suggested 
that a correction of order (3/2) 2 should be included in the dark matter velocity dispersion 
(Turner et al. 1984). This correction is most likely inaccurate since it was computed using 
the global virial theorem, which is not valid for power law profiles that diverge at the origin 
or at infinity. Fitting the observed velocity dispersion as function of radius to tracer mod- 
els of the isothermal sphere profile, it was concluded that ctdm ~ vi os for regions close to 
the center of the galaxies (Kochanek 1993b). This was verified by lensing statistics also, in 
Kochanek (1994). This conclusion is also supported by detailed dynamical modeling of early 
type galaxies (Breimer & Sanders 1993; Franx 1993; Kochanek 1994). There is also a 30% 
uncertainty in the value of F, which we did not take into account. This uncertainty is likely 
to become smaller once more data on galaxies become available. 

In our estimates, effects such as reddening and the luminosity of the lensing galaxy 
(Kochanek 1991) were ignored as there is still not enough quantitative knowledge of their 
results, but it is known that for a faint survey, such as the SDSS, effects due to the luminosity 
of the lensing galaxy must be taken into account. 

It is a pleasure to thank A. A. Sen, and L. Teodoro for the discussions on gravitational 
lensing. O.B. acknowledges the partial support of Fundagao para a Ciencia e a Tecnologia 
(Portugal) under the grant POCTI/1999/FIS/36285. 
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Table 1: SNAP specifications for a two year period of observations. 



Redshift Interval 


z= 0-0.2 z=0.2-1.2 


z=1.2-1.4 


z=1.4-1.7 


Number of SNe 


50 1800 


50 


15 
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Table 2. Parameters for the considered fiducial models. 









A s 


a 


Model I 


0.05 


0.95 


0.83 


0.4 


Model II 


0.05 


0.95 


0.758 


0.0 




Fig. 1. — Expected SNAP confidence regions assuming precise knowledge of M. plus the 
CMBR constraints from Bento et al. (2002b). 




Fig. 2. — Expected SNAP confidence regions for Model I assuming no knowledge of M. 




Fig. 3. — Expected SNAP confidence regions for Model II assuming no knowledge of M.. 
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Fig. 4. — Lensing 68% confidence regions for Model I. The thick lines are the CMBR con- 
straints from Bento et al. (2002b). 
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Fig. 5. — Joint 68% confidence regions for Model I from lensing, SNe and CMBR constraints 
from Bento et al. (2002b). The SNe confidence region is computed assuming an exact M. 
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Fig. 6. — Lensing 68% confidence regions for Model II. The thick lines are the CMBR 
constraints from Bento et al. (2002b). 
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Fig. 7. — Joint 68%CL confidence regions for Model II using both SNe, gravitational lensing 
statistics and CMBR constraints from Bento et al. (2002b). The SNe confidence region is 
computed assuming an exact M. 
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Fig. 8. — Joint 95%CL confidence regions for Model II using both SNe, gravitational lensing 
statistics and CMBR constraints from Bento et al. (2002b). The SNe confidence region is 
computed assuming an exact M. 



